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ABSTRACT 

We have characterized the amount, spatial distribution and kinematics of the 
molecular gas in the merging, double nucleus type 2 quasar SDSS J0025-10 at z =0.30 
using the CO(l-O) transition, based on data obtained with the Australia Telescope 
Compact Array (ATCA). This is one of the scarce examples of quasar host galaxies 
where the CO emission has been resolved spatially at any redshift. We infer a molecular 
gas mass Mh2=(6±1)x10^ M©, which is distributed in two main reservoirs separated 
by ~9 kpc. ~60% of the gas is in the central region, associated with the QSO nucleus 
and/or the intermediate region between the two nuclei. The other 40% is associated 
with the northern tidal tail and is therefore unsettled. 

With its high infrared luminosity Lir=(1.1±0.3)x1012 Lq, SDSS J0025-10 is an 
analogue of local luminous LIRGs and ULIRGs. On the other hand, the clear evidence 
for an ongoing major merger of two gas rich progenitors, the high Ljr dominated by 
a starburst, the massive reservoir of molecular gas with a large fraction still unsettled, 
and the quasar activity are all properties consistent with a transition phase in the 
(U)LIRG-optical QSO evolutionary scenario. We propose that we are observing the 
system during a particular transient phase, prior to more advanced mergers where the 
nuclei have already coalesced. 

We argue that a fraction of the molecular gas reservoir is associated with a tidal 
dwarf galaxy identified in the optical Hubble Space Telescope image at the tip of the 
northern tidal tail. The formation of such structures is predicted by simulations of 
colliding galaxies. 

Key words: galaxiesievolution; galaxies interactions; quasars:individual: SDSS 
J002531.46-104022.2 



1 INTRODUCTION 

CO studies of type 1 quasars (QSOl) at different z have 
provided accumulating evidence that they often contain an 
abundant supply of molecular gas. At redshift z <0.5, the 
range of molecular gas masses spans at least two orders of 
magnitude: some quasars have Mna^lO* Mq significantly 
lower than the Milky Way (few x 10^ M©) and others con- 
tain giant reservoirs >10^'' M©. At the high redshift {z >2), 
quasars contai n Mh , ^several x 10^ to ~10^^ M© (e.g. 
Bertram et al. 12001 Xia et al. |2012| '1. As found for other 



galaxy types, these works show that quasars with the high- 
est infrared luminosities also contain the largest Mhj • 

The investigation of the molecular gas content of type 
2 quasars (QS02jj has been quite active in recent years, 
although almost invariably focusscd on high z >2 objects 
(e.g. Martinez Sansigre et al. 2009 ') and with scarce spatial 
information. 

Regarding studies of the molecular gas content of 
quasars in general, the intermediate redshift range (0.1< 
z <1.5) has remai ned p ractically unexplored until very re- 
cently (Xia et al. |2012| . Krips, Neri & Cox |2012| . Villar- 



* Based on observations carried out with the Australia Telescope 
Compact Array 



^ We refer to QS02 as radio quiet objects, as opposite to narrow 
line radio galaxies 
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Martin et al. 120131 ). This redshift range spans ~60% of the 
age of the Universe, an epoch of declining cosmic star for- 
mation rate (Hopkins & Beacom 2006). Results for QS02 
at redshift z <0.5 h ave be en published only in the la st year 
(Krips, Neri & Cox |2012| . Villar-Martm et al. |2013| ). They 
suggest that QS02 often harbor large reservoirs of molecu- 
lar gas (^severalxlO^ M©), similar to QSOl of similar IR 
luminosity. 

Detailed studies of the spatial distribution of the molec- 
ular gas in quasar host galaxies is of special interest. This gas 
is highly sensitive to the different mechanisms at work dur- 
ing galactic evolution. As such, it retains relic information 
about the global history of the systems. The spatial redistri- 
bution of the molecular gas during the interactions/merger 
processes is determinant on the triggering of the starburst 
and (probably) the q uasar activities (Bessiere et al. l2012l . 
Ramos Almeida et al. l201lh , while powerful negative feed- 
back might be able to counteract this by cleaning and/or 
destroying the molecular gas reservoirs and so quenching 
the star f ormation and quasar activities (Hopkins & Bea- 
com |200i). 

Such studies are very scarce and have also mostly fo- 
cussed on the highest redshift QSOl (e.g. Carilli et al. i2002h . 
The main reasons have been so far the lack in sensitivity 
and the narrow bandwidths of (sub)-mm telescopes as well 
as the major interest on the most distant sources. Spatially 
resolved data on QSO in the local and intermediate redshift 
{z <1 .5) u niverse are also limited in number (e.g. Fe ruglio 
et al. |2010| . Aravena et al. I2OIII , Papad opoulos et al. |2008| . 
Krips et al. I2OO7I . Staguhn et al. l2004l ). The resuhs so far 
indicate that a large fraction of the molecular gas in low 
z quasars (as well as in luminous [10^^ ^Lir/L0 <10^^] 
and ultraluminous [lO^'^ ^Lir/L© <10^''] infrared galaxies, 
LIRGs and ULIRGs) is highly concentrated in the galac- 
tic centers, most frequently in rotating disks and rings of 
<few kpc diameter (e.g. Westmoquette et al. l2012l . Bryant 
& Scoville [1999I) . At intermediate z, some systems have ex- 
tended and massive reservoirs of molecular gas, associated 
wi th com panion objects and / or tid al structures (Aravena et 
al. l201ll . Papadopoulos et aL l2008h . 

We present here a detailed study of the spatial distribu- 
tion of the molecular gas traced by the CO(l-O) transition 
in SDSS J002531.46-104022.2 (SDSS J0025-10 hereafter), a 
QS02 at 2=0.303 selected from Sloan Digital Sky Survey 
(SDSS) database. The CO(l-O) (i/rest=115.27 GHz) transi- 
tion is the least dependent on the excitation conditions of 
the gas, which is crucial for deriving reliable estimates of 
the total molecular gas content, including the wide-spread, 
low-density gas t hat rn ay be sub-thermally excited (e.g. Pa- 
padopoulos et al. l200ll . Carilh et aL boid ). 

SDSS J0025-10 has special interest since it is undergo- 
ing a strong transformation via a major merger event and 
allows to investigate the distribution and dynamics of the 
molecular gas during the process. It moreover shows evi- 
dence for phenomena frequently found in quasars which can 
have a profound impact on the evolution of their hosts: 
mergers/interactions, star formation and nuclear outflows. 
A complete study of this system requires the quantification 
of the molecular gas content and the characterization of its 
spatial distribution and kinematics. This is the purpose of 
this paper. 



1.1 SDSS J0025-10 

Villar-Martin et al. (2011a, 2011b, VMlla and VMUb here- 
after) present a detailed study of this system based on deep 
optical imaging and spectroscopic data obtained with the 
Faint Object FOcal Reducer and low dispersion Spectro- 
graph (F0RS2) on the Very Large Telescope (VLT). Based 



on the high [Oni]A5007 luminosity (log(. 
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=8.73) and 



other criteria related to the optical emission line ra- 
tios and widths , the object was classified as a QS02 by 
IZakamska et al.l (|2003l ). This QS02 is a member of an in- 
teracting system. It has two nuclei located at ^^5 kpc in 
projection. The velocity shift AV=-20±20 km s~^ between 
them implied by the optical spectrum (VMllah . is consistent 
with their relative motion being approximately constrained 
on the plane of the sky and thus, the apparent separation is 
similar to the true separation. One nucleus hosts the quasar 
and the other is forming stars actively. Tidal tails stretch 
to the North and South. Their predominance in the broad 
band optical images shows that they are dominated by stel- 
lar continuum emission, although line emission from ionized 
gas is also detected. 

The optical emission line spectra reveal extended recent 
star formation in the companion nucleus (that sometimes we 
will refer to as nuc2) and the northern tidal tail. Since only 
stars with masses of >10 Mq and lifetimes of <20 Myr con- 
tribute significantly to the integrated ionizing flux, this sets 
an upper limit on the age of the most recent burst of star for- 
mation at these locations. Moreover, compact knots that are 
visible in the tidal tail (see Hubble Space Telescope (HST) 
image in Villar Martin et al. 2012) are reminiscent of star 
clusters and/or tidal dwarf galaxies. The presence of a young 
stellar population (<40 Myr) in the quasar nucleus is also 
conflrmed by fits of the spectral energy distribution (SED) 
of the optical continuum (Bessiere et al. 2013, in prep.). 

The quasar nucleus hosts an ionized outfiow of uncer- 
tain origin (AGN and/or starburst induced), which pro- 
duces very broad emission line components with full width 
half maximum FWHM'-^1300 km s"'^ and blueshifted by '^- 
80 km s ~^ relative to the systemic redshift of the galaxy 
(jVMllbl ). 

We assume Q.a=(I.T, nM=0.3, Ho=71 km s~^ Mpc~^ 
At z =0.303, l"corresponds to 4.45 kpc. 



2 ATCA OBSERVATIONS AND HST ARCHIVE 
DATA 

The observations were performed during 2-7 August 2012 
with the Australia Telescope Compact Array (ATCA), a ra- 
dio interferometer in Narrabri, Australia. The ATCA was 
configured in the most compact hybrid H75 array configu- 
ration. Two 2 GHz bands with 1 MHz channel resolution 
were centered on the redshifted frequency of the CO (1-0) 
line (88.439 GHz), providing redundancy in case technical 
issues would occur with one of the bands. This resulted in 
a velocity coverage of 7000 kms~^, maximum resolution of 
3.5 kms~^ and primary telescope beam (i.e. effective field- 
of-view) of 32 arcsec. Observations were done above an ele- 
vation of 35° and under good weather conditions, with sys- 
tem temperatures ranging between 450 and 800 K (depend- 
ing on the antenna and the elevation of the source) and 
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typical atmospheric seeing fluctuations < 150 /im (phase 
de-co rrelation of < 10%, see iMiddelberg. Sault fc KesteverJ 
uOOa ). The total on-source integration time was 17 hours. 

The phases and bandpass were calibrated every 7.5 min- 
utes with a short {~2 min) scan on the nearby bright cal- 
ibrator PKS 0003-066 (S'88.4GHz = 1.8 Jy at 88.4 GHz), lo- 
cated at 6.4° distance from our target source. Atmospheric 
amplitude variation were calibrated every 30 minutes using 
a paddle scan, and telescope pointing was updated every 
hour, or every time the telescope moved >20° on the sky. 
For absolute flux calibration, Uranus was observed close to 
our target source. 

The off-line data reduction w as don e with the MIRIAD 
software (Sault, Teuben & Wright [l 9951 1. Bad data (includ- 
ing data with internal interference or shadowing of an an- 
tenna, or data taken during weather conditions that in- 
troduced significant phase decorrelation) were discarded. 
PKS 0003-066 was used for both phase as well as time- 
dependent bandpass calibration by interpolati ng the calibra- 
tion solution obtained every 7.5 minutes (see lEmonts et al.l 
I2OIII for details on time-dependent bandpass calibration 
within MIRIAD). Atmospheric opacity variations were cor- 
rected by weighting the data according to their 'above at- 
mosphere' (i.e. paddle-corrected) system temperature. Flux 
calibration was applied using Uranus [ATCA Uranus model 
version Aug 2012], resulting in an absolute flux calibration 
accuracy of 20%. After Fourier transfor mation, we ob tained 
a data cube with robust weighting -|-1 (|Briggslll995l ). The 
data presented in this paper were binned by 10 channels and 
subsequently Hanning smoothed to a velocity resolution of 
68 kms~^, yielding a noise level of 0.7 mjy beam~^ chan~^. 
The synthesized beam-size of the data is 6.30 x 4.39 arcsec^ 
(PA -85.3°). Total intensity images were created by sum- 
ming all the signal (i.e. without setting a noise treshold) 
across the velocity ranges in which CO (1-0) was detected. 
The calculated CO luminosities in this paper have been de- 
rived from these total intensity images (see §3.1). 

No 88.4 GHz radio continuum was detected in the data 
down to a 5(t limit of 0.5 mjy. 

The HST image, which was obtained with the 
ACS/WFC (Advanced Camera for Surveys/Wide Field 
Camera), was retrieved from the HLA (Hubble Legacy 
Archive; programme identification 10880 and principal in- 
vestigator H. Scmitt). The only processing applied was cos- 
mic ray removal. The accuracy of the HST astrometry is 
^0.3" in both coordinates. The accuracy of the CO astrom- 
etry is expected to be sig nificantly smaller (<0.1"; see e.g. 
iPapadopoulos et al.|[200a ). 
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Figure 1. Fit of the near to far infrared SED of SDSS J0025- 
10 based on WISE and IRAS photometry and the ATCA 3mm 
continuum upper limit. Detections and upper limits are marked 
with diamonds and triangles respectively. Different line styles are 
used for the two SEDs discussed in the text. 



band is known to be a complex mixture of stellar and AGN 
related components (scattered and/or direc t AG N light and 
ne bular continuum; e.g. Tadhunter et al. 1201 ll . Vernet et 
al. I2OOII '). To build the SEP we used the SWIRE template 
library (PoUetta et al. 120081 ) which contains 25 templates 
including ellipticals, spirals, starbursts, type 2 and type 1 
AGN and composite starburst -I- AGN. 

The best fit corresponds to the starburst tem- 
plate IRAS 20551-4250 (Fig. 1, dashed line). We infer 
LiR = (1.4±0.4)xlO^^ L0 for this SED. The starburst tem- 
plate NGC 6240, which also produces a reasonable fit, is 
also shown for illustration (Fig. 1, dot-dashed line) because 
this LIRG system is similar in many aspects to SDSS J0025- 
10 (double n uclei, i nteracting local LIRG, Tecza et al. i2000l . 
Tacconi et al. l 19991 ). For this SED, Lir=(8.2±0.2) x lO" Lq. 

With these two SEDs we constrain the infrared lumi- 
nosity LiR= (l.l±0.3)xl0^'^ Lq. Taking the uncertainties 
into account, SDSS J0025-10 is in the transition range of 
LiR values between the LIRG and the ULIRG regimes. 

The fits show that the infrared infrared luminosity of 
SDSS J0025-10 is dominated by a starburst. We calculate 
the star forming rate as SFR = 1.73x 10"^° xLir (M© yr"^) 
— 190±52 Mq yr~^, assuming solar abundances and a 
Salpeter Initial Mass Function (IMF) (Kennicutt 1998). 



3 RESULTS 

3.1 SDSS J0025-10 is a (U)LIRG 

We have constrained the infrared (IR) luminosity, Lir (in- 
ferred from the 8-1000 /im wavelength range) by fitting the 
source spectral energy distribution (SED) defined by the 
Wide-Field Infrared Explorer (WISE; 3.3, 4.6, 11.6, 22.1 
yitm) and Infrared Astronomical Satellite (IRAS; 60 and 100 
/im) photometric measurements, as well as the upper limit 
on the ~3 mm continuum measured with the ATCA data 
(see Fig.l). Optical photometry has not been used, since this 



3.2 Mass and spatial distribution of the molecular 
gas 

We measure a total CO (1-0) luminosity in SDSS J0025-10 
of L'co=(7.5±1.5)xl0® K km s"^ pc^. Lqo is calculated as: 



L'co =3.25 X 10^(1^^) ( 



Dl , 2 / frost N - 2 / , . N - 1 

) (fttjz) (1 + 2) 



Jy km/s Mpc GHz 

where Ico = Sco AV is the integrated CO (1-0) line 
intensity in Jy km s~^, Dl is the luminosity distance in 
Mpc and :^rest = 115.27 GHz, is the rest frame frequ ency o f 
the CO (1-0) transition (Solomon & Vanden Bout 120051 '). 
With this L'co, the system follows the L'co vs. Lir (or 
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Figure 2. l-dimcnsional spectrum of SDSS J0025-10, extracted 
against the QS02 nucleus (2.3" x2.3"). Notice the double horned 
profile, with Gaussian fits visualizing the two components de- 
scribed in the text. The bottom plot shows the residuals after 
subtracting the two Gaussian profiles from the spectrum. Veloc- 
ity values are relative to the peak of the 'blue' component, which 
we assume as the systemic velocity (z=0.3031±0.0001, see text). 



Lfir) correlation define d by difTerent types of galaxies 
(Solomon & Va nden B out l2005l ). includin g QS 1 and QS02 
(Bertr am et al. |2007| . Villar-Martin et al. |2013| . Krips, Neri 
& Cox l2012l ). The molecular gas mass is then calculated as 
Mh2 = axL'co = (6±l)xlO^ Mq, assuming a=0.8 Lq (K 
km s~^ pc^)~^. This value of a is often used in studies 
of (U )LIRGs and active galaxies (e.g. Downes & Solomon 
1 19981 ). although a range a ^0.3-1.3 is possibl e (e.g . Solomon 
& van den Bout |2005| . Sanders & Mirabel [1999 ). Similar 
amounts of molecular gas have been measured for quasars 
(both type 1 and type 2) of simil ar infr ared luminosity and 
redsh ift (e.g. Krips, Neri & Cox |2012| . ViUar-Marti'n et al. 
l2013l ). 

The CO (1-0) line shows a double horned profile (Fig. 2) 
with two kinematic components. The blue component has 
z=G. 3031±0.0001 (which we assume as the systemic redshift 
Zsys) and FWHM=140±25 km s"^ The red, fainter compo- 
nent is shifted by -f 160±10 km s"^ and has FWHM=80±40 
km s~^. A negative depression (around v = 500kms~^) is 
also visible in Fig. 2. This feature appears at -2.8a (when 
integrated over its full extent). Fig. 2 (bottom) shows that 
it may be merely a prominent noise peak in the data. Since 
this negative depression is not noticeably caused by ex- 
ternal effects (e.g. weather or phase calibration), it may 
also indicate that there are low-level systematics affect- 
ing the A TCA 3mm (similar t o low-level systematics de- 
scribed by lEmonts et al.l l201ll that occur at the extreme 
edge of the A TCA 7mm band, or possible baseline artifacts 
mentioned bv lNorris et al.l (|2013l ) in early wide-band 3mm 
ATCA data). Regardless, the integrated blue and red part 
of the double horned CO profile are detected at the IOct and 



7(T significance level respectively and thus unambiguously 
stand out above the noise, indicating the reliability of the 
CO emission discussed in this paper. 

The spatial distribution of CO (1-0) in different chan- 
nel maps is shown in Fig. 3 overlaid on the HST optical 
image. The CO(l-O) emission detected at ^3a level spans 
a range of velocities ~[-110, -1-230] km s"'^ relative to Zsys- 
The CO (1-0) is spatially extended and accumulated in two 
main reservoirs separated by ~2"or ~9 kpc (blue and red 
contours in Fig. 3, first panel) which we will name J0025A 
and J0025B. The similarity within the errors of the velocity 
shift between them, compared to the two kinematic compo- 
nents isolated on the 1-dimensional spectrum suggests that 
J0025A and J0025B are responsible for the double horned 
profile discussed above. J0025A and J0025B have CO lumi- 
nosities of (4.3±0.9)xl0^ and (3.1±0.9)xl0^ K km s~^pc^ 
respectiveljo, which correspond to Mh2 =(3.4±0.7) x 10^ and 
(2.5±0.7)xl0^ Mq. 

In Fig. 4 we show the location of the peaks of the CO(l- 
0) emission in J0025A and J0025B overlaid onto the HST 
image. These locations have been derived by fitting a Gaus- 
sian profile to the integrated CO (1-0) emission. Within the 
uncertainty, this Gaussian profile could be represented by 
a point source with the size of the synthesized beam. The 
locations of the CO (1-0) peaks in the different channel maps 
of Fig. 3 are also plotted, including the uncertainty due to 
fitting and astrometric errors (see caption for details) . From 
Fig. 4 it is clear that the CO(l-O) emission in J0025B is as- 
sociated with the northern tidal tail while J0025A is prefer- 
entially associated with the QS02 nucleus and/or the region 
between both nuclei. 

We do not detect any CO (1-0) emission associated 
with the companion nucleus nuc2. However, the star form- 
ing activity it habours (§1.1) requires a loc al supp ly of 
molecular gas. The H/3 luminosity from nuc2 (|VMlla|) im- 
plies a line reddening corrected SFR^8 M© yr~^, assuming 
SFR(M0 yr"^) = 2.3x10"*^ - ^^ for solar abundances 
and a Salpeter IMF (Kennicutt ll998l )FI This in turns im- 
plies LiR,i„c2~4.6xlO^'' Lq, which would be a lower limit if 
part of the st ar form ation were totally obscured (Rodriguez 
Zauri'n et al. ^^j. This is -8% of the total Lir of the 
system. If nuc2 follows the Lir vs. L'co correlation defined 
by st ar for ming galaxies (see Fig. 8 in Solomon & van den 
Bout|2005|), it roughly implies that it contains >8% of the 
total molecular gas or >5xl0* Mq. The 3a upper limit for 
a CO (1-0) non-detection with FWHM = 140 km s"^ in our 
dataisL'co= 1.0x10^ K km s"^pc^ or Mh^ = 8x10** Mq. 
Although the calculations are highly uncertain due to the 
large scatter of the Lir vs L'co correlation and the uncer- 
tainty on the amount of obscured star formation in nuc2, 
they imply that its CO(l-O) emission likely falls below the 
detection limit of our data (or otherwise a low-level signal 



Values have been derived from the total intensity image in 
Fig. 3, but were corrected for line contamination and missed low- 
level flux based on the Gaussian fitting of the profiles (Fig. 2). 
This did not significantly change the values for the blue compo- 
nent, but increased both the estimated value and the associated 
uncertainty of the fiux of the red component by 10%. The errors 
also reflects the 20% uncertainty in absolute flux calibration (§2). 
^ For comparison, the SFR over the entire Milky Way is several 
Mq yr-i (e. g. Robitaille & Whitnev [JoToh 
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Figure 3. The total intensity image of both the blue (J0025A, -90 s£ v < 79kms-l) and red (J0025B, 79 s£ v C 249 kms"!) part of 
the CO(l-O) signal is shown in the first frame, followed by the channel maps of the CO emission. Contour levels of the total intensity 
map: 0.27, 0.41, 0.54, 0.68, 0.81 Jy bm"! X kms~i; contour levels channel maps: -3.0 (dashed), 3.0, 4.5, 6.0, 7.5, 9.0 (solid) X a, with 
(T=0.7 mjy bm~^. Velocity values are relative to Zsys = 0.3031 it 0.0001. The location of the QS02, companion star forming nucleus and 
the northern tidal tail are indicated in the second panel. A color version of this figure can be seen in the electronic form. 



would likely be hidden within the large ATCA beam by the 
dominant emission from the quasar nucleus and/or interme- 
diate region). New observations with higher spatial resolu- 
tion and sensitivity are required to both determine the exact 
location of the CO in J0025A and measure (or set reliable 
limits on) the CO content of the star forming companion 
nucleus. 



4 DISCUSSION 

4.1 SDSS J0025-10 vs. local (U)LIRGs 

SDSS J0025-10 is a (U)LIRG double nuclei merging sys- 
tem with all observed properties consistent with those of 
numerous local [z <0.05) luminous LIRGs (severalxlO^'^ < 
LiR<10^'^ Lq) and ULIRGs. It shows clear morphological ev- 
idence for an ongoing major merger: the double nucleus, the 
tidal tails, the inner spiral structure (which surrounds the 
QSO nucleus clockwise and apparently becomes the north- 
ern tidal tail) and the stellar clusters and/or tidal dwarf 



galaxies observed along this structure (Fig. 4) are all mor- 
phological features o ften o bserved in ULIRGs and lu minous 
LIRGs (Haan et al l201ll . Miralles Caballero et al. 1201 ll , 
Monreal Ibero et al. 120071 ). As in these systems, the star for- 
mation and probably the AGN activ ity are triggered by the 
interaction (see Al onso-H errero l2013l for a recent review; see 
also Bessiere et al. 120121 ) . 

The IR luminosity, Lir=(1.1±0.3)x10^^ Lq, which is 
dominated by the starburst component, the large molec- 



ular gas content Mh2=(6±1)x10" M 

Lo M 



01 



=190±52 



the SFR= 

Mo yr"\ the star formation efficiency SFE=-^IB- = 183±62 

ivi0 ^ and the gas exhaustion time scale tsf = g^^ = 
32±11 Myr are all in the range of lum inous LIRGs and 
ULIRGs (Solomon & van d en Bout l2005l . Sanders & Mirabel 
1996, Alonso-Herrero l2013l ). 

In SDSS J0025-10 one nucleus hosts the AGN and both 
harbour star formation. About 20-30% of local LIRGs have 
an optically identified AGN, with this percentage increasing 
with the LiR and becoming ~70% for ULIRGs (e.g. Yuan, 
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Kewley fc S anders [JOlol . Nardini et al. |2010| . Veilleux et al. 
1 19991 . 119951 ). On the other hand, a large fraction (~63 %) 
of local lu minous LIRGs and ULIRGs have double nuclei 
(Haan et al. 120111 ) with a variety of scenarios. In some cases 
both nuclei host star formation and AGN activities (being 
NGC6240 the most famous example; e.g. Tecza et al. 2000, 
Tacconi et al. 1999); in others star formation is detected in 
both nuclei, but none seem to harbour AGN activity or show 
a composite AGN+starburst spectrum, like SDSS J0025-10 
(e.g. Yuan, Kewley & Sanders [JOlol ). 

The CO gas in SDSS J0025-10 is distributed in two 
main spatial components, J0025A and J0025B of masses 
(3.4±0.7)xl0'' and (2.5±0.7)xl0^ Mq respectively. They 
are separated by 2"or ~9 kpc (projected distance) in space 
and 160 km s~^ in velocity. J0025A is preferentially associ- 
ated with the quasar nucleus and/or the intermediate region 
between the two interacting nuclei, while J0025B is associ- 
ated with the northern tidal tail. CO (1-0) emission from the 
companion nucleus is not unambiguously detected, although 
the active star formation it hosts implies that there must be 
a local reservoir of Muai^severalxlO* M©. The CO gas in 
local (U) LIRGs is usually highly concentrated in the central 
region towards the core of the merg er (r <1 kpc, e.g. Bryant 
& Scoville [l999l . Sanders & Mirabel [l99^ ). In addition, a sig- 
nificant fraction of the gas is sometimes spread over scales 
of 10 kpc or more. This is possibly the case of SDSS J0025- 
10: J0025A might be the compact, central concentration (al- 
though higher spatial resolution observations are required to 
constrain its exact location and compactness), while J0025B 
is the more extended gaseous component. 



4.2 SDSS J0025-10 as a (U)LIRG-optical QSO 
transition object 

Both theory and observations support that the merger of 
two gas rich progenitor galaxies form a ne w, mo re massive 
elliptical galaxy (e.g. Toomre & Toomre Il972l . Barnes & 
Hernquist [l99y) . As a consequence of the gravitational in- 
teraction, gaseous dissipation funnels gas into the the grav- 
itational center and triggers intense star formation, which 
is manifested in the (U)LIRG phenomenon. The gas inflow 
could also fuel the formation of an AGN (e.g. Hopkins et al. 
120081 . 120061 ). The discovery of massive reservoirs of molecular 
gas shifted far away from the quasar location has prompted 
some authors to suggest that gas poor/gas ric h gala xy merg- 
ers c an also trigger quasar activ ity (R iechers l2013l . Aravena 
et al. l201ll . Papadopoulos et al. l2008l ). 

We are observing SDSS J0025-10 when the galaxies are 
deeply interpenetrating and approaching the final coales- 
cence (r <a few kpc). Simulations of gas dynamics and star- 
bursts in major gas rich mergers with or without quasar 
activity (e.g. Mihos & Hernquist 1996, Hopkins et al. |2005| ) 
show that at this stage the morphologies show distorted, 
irregular isophotes and extended tidal tails, as observed in 
SDSS J0025-10. The models also predict that at this moment 
the gas experiences the strongest levels of inflow so that large 
amounts of gas accumulate in the central region and very in- 
tense star formation activity is triggered as a consequence. 
The duration of the starburs t is ex pected to be very short 
~50 Myr (Mihos & Hernquist 1 1994 ) . The presence of such a 
young stellar population is confirmed in SDSS J0025-10 (see 
§1-1)- 



ISanders et al.l l|l988al ) proposed that there is a natu- 
ral evolution from ULIRG to QSO: ULIRGs are the initial 
dust-enshrouded stage of quasars, which will become op- 
tical quasars once the nuclei shed the obscuring d ust vi a 
feedback mechanisms (see also Hopkins & Beacom 120051 ). 
Some authors have expanded the ULIRG-QSO evolution- 
ary scenario accommodating luminous LIRGs as a phase 
that appea rs earl ier in mergers than the ULIRG phase (e.g. 
Haan et al. l201ll . Yuan, Kewley & Sanders [2010). A review 
on the observational and theoretical studies, with pros and 
cons of these sc e nario s can be found in the recent paper by 
iRothberg et al.l l|2013l ). Although the merger driven evolu- 
tionary sequence (U)LIRG-QSO probably is too simplistic 
(Veilleux et al. .2002 ), both theory and observations support 
that at least a signiflcant fraction of optical quasars might 
evolve from a prior (U)LIRG phase. 

In such scenario, SDSS J0025-10 is a strong candidate to 
be a (U)LIRG-optical QSO transition object. All the prop- 
erties investigated here are expected for such a transition 
phase: high infrared luminosity powered by intense star for- 
mation, a rich molecular gas reservoir, strong evidence for a 
major merger of two gas rich progenitors and optical id en- 
tiflcation of quasar activity (Sanders et al. Il988al . ll988bl ). 

The nuclear separation in SDSS J0025-10 (~5 kpc) is 
consistent with typical values of double nuclei LIRGs and 
signiflcantl y larg er than in more advanc ed m ergers (e.g. 
Haan et al. l201ll . Yuan, Kewley & Sanders [2010|), which are 
more frequently found in ULIRGs. This and the fact that 
a large fraction of the molecular gas is still unsettled show 
that the system is in an transient stage prior to the final 
coalescence of the nuclei. 

The identificati on of an ionized outfiow associated with 
the quasar nucleus (|VMllbl ). prompts the interesting possi- 
bility of tracing its signature in both the neutral and molec- 
ular phases to fully characterize its mass, geometry and en- 
ergetics. This information would allow to evaluate whether 
the outflow is powerful enough to clear up completely the 
material surrounding the quasar nucleus, quench the star 
formation activity and convert SDSS J0025-10 into an op- 
tical quasar before becoming a non active m assive elliptical 
galaxy (e.g. di Matteo, Springel & Hernquist |2005| ) . 



4.3 A large reservoir of molecular gas at the tip 
of the tidal tail 

Merger simulations of gas rich galaxies show that, besides 
the infalling gas which concentrates in a compact central 
region of less than 1 kpc in size, part of the gas survives the 
merger (Barne s fc H ernquist 1996; Mihos & Hernquist 1996; 
Hopkins et al. 120081 ) . This is primarily material which has 
been moved to large radii temporarily, either blown out by 
a combination of supernova and AGN feedback or thrown 
out in tidal tails. The lack of spatial information prevents us 
from making a confldent elucidation of the nature of J0025A 
and J0025B but given their spatial location, it is reasonable 
to propose that J0025A is (or will become) the compact 
component while J0025B, associated with the northern tidal 
tail, is "surviving" material. 

In Fig. 4, it seems that the locations of the peak of 
the CO (1-0) emission in the various channel maps traces 
the tidal tail, from the inner spiral-like structure between 
the nuclei to the tip of the tail, maybe beyond. This sug- 
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gests that there is molecular gas along the tail. According 
to the model predictions, this surviving gas will then slowly 
settle back into a small, rotationally suppo rted e mbedded 
disk (Hopkins et al. 2008, Bournaud & Due 120061 ') ~a cou- 
ple X 10* yr after the coalescence. This gas is not expected 
to contribute much to the galaxy stellar mass (~a few per 
cent), even for the most gas rich mergers. Thus, although 
a large reservoir of gas is still falling towards the center, it 
seems unlikely that it will produce a significant enhancement 
on the build up of the galaxy stellar mass. 

On the other hand, the centroid of J0025B is coincident 
with a prominent feature of elliptical shape located at the tip 
of the northern tidal tail and clearly identified in the HST 
image (Fig. 4). A CO accumulation at a tidal tail's tip is a 
promisi ng ti dal dwarf galaxy (TDG) candidate (Bournaud 
& Due 120061 ). A TDG is a self-gravitating entity that was 
made up from the debris of a galaxy interaction. 

Indeed, the properties of the optical feature strongly 
support this scenario. The optical size along the major 
axis D =1.4±0.1 kpc, the effective radius Reff =0.70±0.06 
kpc, the observed H/3 luminosity L(H/3)=3.6xlO*° erg s~^ 
IVMllaf ) are all consiste nt wit h values measure d for TDGs 
(e.g. Monreal Ibero et al. l2007l . Bournaud & Duc l2006l '). The 
J0025B CO (1-0) line luminosity is >50 times higher (and 
therefore implies >50 times more molecular gas mass) than 
typic ally measured for TDGs (e.g. Braine, Due & Linsfeld 
I2OOII ). but this most probably reflects the CO (1-0) contam- 
ination of J0025B by emission from other regions. 

The dynamical mass o f the TDG can be estimated as 
(Monreal Ibero et al. l2007f) : 



Mdyn(M0) = constant x 10 — ^ ( -r) 

kpc km s 

The constant value ranges between 1.4 and 2.2. Follow- 
ing Monreal Ibero et al. (|2007|) we adopt a value of 2.09. 
This calculation is highly uncertain due to the uncertainty 
on the exact mass distribution of the object and specially 
on the velocity dispersion (t=FWHM/2.35. In spite of this, 
an estimation of the ranges allowed can be useful to sup- 
port or discard the TDG scenario. The FWHMco=165±85 
km s~^of J00255j is not a reliable tracer of the dynamical 
mass, since, as we have seen, it is contaminated by emission 
form other regions. On the other hand, the optical emission 
lines, which are much less affected by this problem, have 
FWHM(H/3)s;i40 km s-\ FWHM([OIII]A5007)=330±10 
kms~^ (VMlla ) and FWHM([OI]A6300)=140±20 km s^V 
Such discrepancy for [OIII] might indicate that it is broad- 
ened by some turbulent mechanism (e.g. winds). It is reason- 
able to assume FWHM'~140 km s~^or a ~60 km s"'^ (notice 
that this value also locates the system with other TDGs in 
the Reft vs. a plane; e.g. Monreal Ibero et al. 120071 ). In such 
case, Mdyn ~5xl0^ Mq. Dynamical masses as large as this 
have been measured fo r othe r TDG (Monreal Ibero et al. 
I2OO7I . Bournaud & Due 120061 ). At the lowest extreme, even 
unexpectedly small a values (e.g. 10 km s~^) would result 
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^ This value is larger than 80±40 km s"'^ quoted in §3.2 for 
the red fainter component of the double horned CO profile. The 
reason is that the total CO emission associated with J0025B has 
been considered, rather than only the inner 2.3" x2.3" pixel used 
to build the 1-dimensional spectrum discussed in that section. 



0''35'"31.6'' 31.5" 31.4'= 31.3'= 

Right Ascension (J3000) 

Figure 4. Overview of the location of the CO(l-O) peak emission 
against features seen in the HST image (grey-scale). The large '-I-' 
signs indicate the location of the peak of the CO(l-O) emission in 
J0025A (blue) and J0025B (red) derived from the total intensity 
image of Fig. 3 (see text for details). The region covered by the 
size of the '+' symbol indicates the uncertainty in the exact loca- 
tion as a result of uncertainty in the fitting of the CO(l-O) peak 
combined with astrometric uncertainty (§2). The crosses indicate 
the location of the peak of the CO(l-O) emission in the various 
channel maps of Fig. 3 (dark blue — >■ dark red for lowest — > high- 
est velocities). The size of the cross scales with the significance 
level at which the CO(l-O) is detected, while the dashed ellipses 
indicate the uncertainty in the exact location (from the CO(l-O) 
fitting combined with astrometric errors) . A color version of this 
figure can be seen in the electronic form. 



in Mdyn >10* M©, stiU typical of TDGs. Thus, in spite of 
the uncertainties, we confirm that the dynamical mass of 
the optical tidal feature associated with J0025B centroid is 
in the range of TDGs. 

The formation of such massive self-gravitating conden- 
sations of matter in the outer regions of tidal tails is allowed 
by the simulations, without strong constraints on th e galac- 
tic encounter parameters (Bournaud & Due |200J). These 
condensations can survive longer (>1 Gyr) than the inner 
regions of the tails. According to these simulations, only the 
objects formed at the tip of tidal tails can become long- 
lived, massive, dwarf satellite galaxies. Whether this will be 
the case depends on its stability against its own internal 
motions and against tidal forces from the parent galaxies. 



5 SUMMARY AND CONCLUSIONS 

We have characterized the amount, spatial distribution and 
kinematics of the molecular gas in the optically selected type 
2 quasar SDSS J0025-10 at z =0.30, which is a merging sys- 
tem with a double composite nucleus (AGN-I-HII). For this, 
we have used the CO (1-0) transition based on data obtained 
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with the Austraha Telescope Compact Array (ATCA). This 
is one of the scarce examples of quasar host galax;ies where 
the CO emission has been resolved spatially at any redshift. 

The system contains Mh2=(6±1)x10® M©, which is 
distributed in two main reservoirs separated by ~9 kpc. 
~60% of the gas is in the central region, preferentially asso- 
ciated with the QSO nucleus and/or the intermediate region 
between the two nuclei. The other 40% is associated with the 
northern tidal tail and is therefore unsettled. 

Based on WISE and IRAS photometry, we constrain 
the infrared luminosity Lir=(1.1±0.3)x10^^ Lq and show 
that it is dominated by a starburst, rather than the AGN. 
SDSS J0025-10 is an analogue of local [z <0.05) luminous 
LIRGs and ULIRGs. The clear morphological evidence for 
an ongoing major merger, the large molecular gas content, 
the star forming rate SFR=190±52 L© M0"\ the star 
formation efficiency SFE=183±62 L0 M©"^ and the gas 
exhaustion time scale rsp = 32±11 Myr are all in the range 
of these local systems. 

SDSS J0025-10 is the result of an ongoing major merger 
of two gas rich progenitors. This, together with the high 
infrared luminosity dominated by a starburst, the mas- 
sive reservoir of molecular gas and the quasar activity are 
all properties consistent with a transition phase in the 
(U)LIRG-optical QSO evolutionary scenario. The relatively 
large nuclear separation (~5 kpc) and the existence of large 
amounts of unsettled molecular gas (probably surviving gas 
from the merger), moreover suggest that we are observing 
the system during a particular transient phase, prior to more 
advanced mergers where the nuclei have already coalesced. 

We propose that at least part of the extended, unsettled 
reservoir of molecular gas is associated with a tidal dwarf 
galaxy at the tip of the northern tidal tail identified in the 
HST image. The formation of tidal dwarf galaxies in the 
outer regions of tidal tails is predicted by simulations of 
colliding galaxies. 
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